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Abstract 

Recent advances in computing subsonic flow 
have been applied to helicopter configurations 
with various degrees of success. This paper is a 
comparison of two specific methods applied to a 
particularly challenging regime of helicopter 
flight, very low speeds, where the interaction of 
the rotor wake and the fuselage are most 
significant. Comparisons are made between 
different methods of predicting the interactional 
aerodynamics associated with a simple generic 
helicopter configuration. These comparisons are 
made using fuselage pressure data from a Mach- 
scaled powered model helicopter with a rotor 
diameter of approximately 10 feet. The data 
shown are for an advance ratio of 0.05 with a 
thrust coefficient of 0.0066. The results of this 
comparison show that in this type of complex 
flow both analytical techniques have regions 
where they are more accurate in matching the 
experimental data. 

Introduction 

Rotorcraft configurations have always presented 
a challenge to the accurate prediction of vehicle 
aerodynamic performance. Complex lifting 
surfaces operating in a characteristically 
unsteady environment present challenge 
enough, but coupling this flow with the shapes 
common to helicopter fuselages amplifies the 
complexity. Often the prediction of isolated rotor 
aerodynamics is coupled using superposition 
with linear aerodynamic fuselage analyses or 
measured isolated fuselage data. Previous 
studies (reference 1 and 2) of complete 
helicopter configurations have shown 
weaknesses in the linear superposition 
assumptions commonly used in the design 
cycle. However, accurate models for the 
complex aerodynamic interactions between the 
rotor and the fuselage have not been developed 
as general tools available to the helicopter 
designer. 

Non-linear interaction effects arise in the 
aerodynamics of helicopter configurations in 
several cases. Among these cases clearly the 
wake of the rotor affects the fuselage onset flow. 
The wake geometry in most inflow models is 
assumed to be undisturbed by the fuselage. 
The wake does distort due to the presence of 
the fuselage. This distortion increases as the 
wake skew angle decreases (at lower speeds 


where the wake washes over the body). The 
influence of the fuselage on the inflow to the 
rotor is also potentially significant. The additional 
inflow distortion to the presence of the fuselage 
produces a change in the aerodynamics of the 
blades. This effect changes the strength of the 
shed wake and contributes to additional 
distortion of the wake. 

Notation 

Ct Thrust Coefficient, thrust/(density x disk 
area x tip speed xtip speed), where thrust 
is the force normal to the plane of the rotor 

Cp Pressure Coefficient, (pressure - free 

stream pressure)/free stream dynamic 
pressure 

R Radius of the rotor, 62 in (1 .574 m ) 

X Downstream length, in. 

Y Lateral distance, in. 

Z Vertical length, in. 

A1 s Longitudinal Flapping angle, degrees 

Experimental Data 

The Langley 14- by 22-Foot Subsonic 
Tunnel 

The Langley 14- by 22-Foot Subsonic Tunnel is 
a closed-circuit, single-return, atmospheric wind 
tunnel (figure 1). In 1970 the unusual test 
requirements associated with V/STOL and 
rotorcraft aerodynamic research led to design 
and construction of this tunnel. The tunnel has a 
test section that can be operated in a variety of 
configurations: closed, slotted, partially open, 
and open. The closed test section is 14.5 ft high 
by 21 .75 ft wide by 50 ft long with a maximum 
speed of about 338 feet per second (fps). The 
open test section configuration, which has a 
maximum speed of about 270 fps, is formed by 
raising the ceiling and walls to form a floor-only 
configuration. The tunnel may be configured 
with a moving-belt ground plane and a floor 
boundary-layer removal system at the entrance 
to the test section for ground-effects testing. 
During this investigation the tunnel was 
configured with the walls and ceiling in the raised 


1 



position to improve the quality of the low speed 
flow. 



Figure 1 : Langley 14- by 22-Foot Subsonic 
Tunnel 



Figure 2: GRMS and fuselage in Langley 14- by 
22-Foot Subsonic Tunnel 


The rotor system tested consisted of four 
rectangular blades on an articulated hub. The 
blades have a linear twist of -8.0 degrees from 
the center of rotation to the tip. The chord of 
4.25 inches and radius of 62.0 inches gives the 
system a solidity of 0.087. 


The tunnel is equipped with an on-line static data 
reduction system that can display computed 
aerodynamic coefficients with interactions and 
wall interference corrections in real time. The 
tunnel has support, drive, and instrumentation to 
facilitate powered rotorcraft testing and has been 
used for rotorcraft investigations since its 
inception. In addition, the tunnel has flow- 
visualization and acoustic testing capabilities. 

Rotor Test System 

The rotor test system used for the experimental 
data reported here is built on a generic test 
system developed at the 14- by 22-Foot 
Subsonic Tunnel. This test system, the General 
Rotor Model System (GRMS), consists of two 
synchronous electric motors, a combining 
gearbox, collective and cyclic blade pitch 
controls and a four-bladed articulated hub 
mounted on a six-component strain-gauge 
balance. It also includes a fuselage skin mounted 
on a separate, similar balance. These two six- 
component strain-gauge balances are used to 
provide independent measurement of the rotor 
and fuselage aerodynamic loads. The system as 
tested is shown in figure 2 installed in the wind 
tunnel. 


The shape of the fuselage is designed to be 
representative of a wide range of helicopter 
fuselages without being specific to any one. The 
fuselage can be described by a set of super- 
ellipse equations that simplifies development of 
computer models. The geometry of this fuselage 
is described in references 2 and 3 and is shown 
in figure 3. 
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Figure 3: ROtor Body Interaction (ROBIN) 


fuselage (units of nominal rotor radius) 

Test Procedure 

Fuselage surface pressure measurements were 
made at several flight conditions (reference 3). 
Although these conditions ranged in lift, forward 
speed, and propulsive task this comparison will 
be limited to a single representative low-speed 
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condition. The advance ratio for this 
representative condition was 0.05. A non- 
dimensional thrust coefficient of 0.0066 was set 
with a shaft angle of -2.0 degrees relative to the 
free-stream. The rotor tip speed was 649 feet 
per second. The rotor cyclic controls were 
trimmed to reduce once-per-revolution flapping 
of the rotor blades with respect to the rotor shaft 
to -0.8 degrees of longitudinal flapping and 0 07 
degrees of lateral flapping. 

Pressure measurements were made using an 
array of scanning pressure transducers. 

Pressure taps were arranged in strips of taps at 
constant X stations. Only the average pressure 
values could be obtained with the 

instrumentation available during this 
investigation. The helicopter model is mounted 
on a sting support system that fits to the bottom 
of the fuselage and trails the model. 

Computational Methods 

Two methods of computing the interaction of the 
helicopter rotor and its wake with the fuselage 
will be shown in this work. The first method is 
based on techniques that model the flows 
associated with the helicopter. This method has 
been developed by Mil using models developed 
for the analysis of helicopter flows (reference 4). 
The second method shown is based on the 
solution of the Navier-Stokes equations for the 
flow. The effect of the rotor in this method is 
approximated using a “pressure jump” boundary 
condition at the rotor disk. This method has 
been developed by the Aeroflightdynamics 
Directorate (AFDD) of the US Army (reference 5). 

Mil Method 

This work presents some results of the study of 
flow over a helicopter fuselage at the low 
horizontal velocity condition with the influence of 
the main rotor . In these calculations allowances 
were made for wake contraction and for variability 
of convection velocity of free vortices. An 
assumption was made that the free vortices were 
moving down with a velocity averaged over 
azimuth. However this average velocity varies 
with blade radius and according to the vertical 
distance from the rotation plane of the rotor. 
Other details of this method are described in 
reference 4. 


The method shown here is computed in 
sequence on personal-computer class of 
workstations. As such, this method could be 
incorporated in a comprehensive helicopter 
analysis to be a computationally efficient method 
for simulation of the aerodynamic influence of a 
rotor-fuselage configuration. Solutions from this 
method have been provided without an estimate 
of the actual computer time required. Six main 
sub-programs are used in this method. 

1 . The main program that starts and links to the 
rest of the programs. 

2. A sub-program to provide the discrete 
geometry of the fuselage surface and its wake in 
the case where the flow separates at the aft end 
of the fuselage. 

3. A sub-program for the approximation of the 
wake shape when separating at the aft of the 
fuselage using the method described in 
reference 6. 

4. A sub-program for calculating the flow over the 
fuselage. 

5. A sub-program to determine the velocities 
induced by the fuselage and its wake at any 
point in space. 

6. A sub-program to calculate the streamlines on 
the fuselage surface. 

The fuselage geometry is broken down into 
panels using an automated procedure. The 
resulting panels for the Mil calculation are shown 
in figure 4 along with the experimental presure 
measurement stations. Interpolation is required 
between the computed surface pressures at the 
panel centers to the locations of the measured 
pressures. 
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Figure 4: Panels used for Mil fuselage model 

First, an isolated rotor was considered to choose 
values of collective pitch and longitudinal 
deviation of the swashplate that produce the 
necessary thrust coefficient, Ct, and the flapping 
angle, Als. Then the velocity induced by the 
rotor on the fuselage was calculated and 
pressure distribution was defined. 

The calculations of the flow over the fuselage 
were carried out in the velocity field induced by 
the rotor and averaged over time. To determine 
the time-averaged influence of the main rotor, 
the entire rotor wake was computed at three 
azimuths with an interval of 30 degrees. Since 
there are four blades in this rotor system, this 
average is equivalent to a full rotation. 

At the next stage the effect of the fuselage upon 
the rotor was defined. The collective pitch at the 
rotor and the longitudinal deviation of the 
swashplate were slightly corrected, and the 
pressure distribution over the fuselage was 
redefined. It is not necessary to continue the 
iterative process due to the small influence of 
the fuselage on the rotor at the given condition. 
The calculation was then stopped. 

AFDD Method 

The choice of a Navier-Stokes code for this work 
was made to insure that regions of the fuselage 
that would naturally separate could be predicted 
without a priori knowledge of the geometry of 
this separation. An incompressible Navier- 
Stokes code (reference 7) is the basis for the 
method used in this study since the relative 
speeds seen by the fuselage are 
incompressible. The use of a Navier-Stokes 
code for rotorcraft application requires either a 
detailed grid system that describes the rotor 
system and relative motion of the rotor blades or 
a model for the lifting rotor. For the work 


described here, the rotor system was modeled 
by a pressure discontinuity at the rotor disk. The 
compressible effects of the rotor blade are 
accounted for using the blade element theory 
and tabular values for the lift and drag of the 
airfoil at each rotor section. 

The pressure discontinuity at the rotor disk was 
computed iteratively with the solution to the flow 
field. This pressure value was computed using 
blade element theory and the current values of 
the flow field at each location on the rotor disk. 
This computation also included blade pitch trim 
to attain the desired rotor thrust and moment 
balance. 



Figure 5: AFDD Grid for ROBIN and Rotor 


The grid system that was used for the calculation 
of the helicopter flow field consists of several 
overset grids. Interpolation of the flow field 
properties in the overlapping regions of the grids 
was used as the boundary condition for the inner 
grids. A global outer grid was used to prescribe 
the actual flow condition to the code. The 
overlapping inner grids are shown in figure 5. 

The relative cost of doing a complete Navier- 
Stokes simulation of helicopter configurations at 
low speed is high. Over 60 Cray C-90 hours were 
required for this case. 

Results and Discussion 

The character of the aerodynamics at low speeds 
can be seen in the flow solution obtained using 
the AFDD method. Figure 6 shows computed 
flow at the flight condition studied. Since the 
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AFDD model uses steady flow, streamlines 
(shown in sub-figure 6(a)) can be computed to 
show the general directions taken by the flow. 
Shading on the fuselage in sub-figure (a) 
indicates changes in surface pressure. The 
stagnation of the rotor wake over the nose of the 
fuselage is clearly shown in this figure. In sub- 
figure 6(b) the concentration of vorticity in the 
flow is shown by means of the component of 
vorticity in the X-Z plane. This figure shows the 
effective envelope of the rotor wake and the 
influence of the fuselage on this shape. The 
wake skew angle computed by momentum for 
this condition is approximately 45 degrees. From 
the figure, the leading edge of the wake 
envelope is approximately 50 degrees with 
respect to the rotor disk normal. The trailing 
edge of the wake envelope (although clearly not 
fully developed) is less than the leading edge, 
approximately 42 degrees. Beyond the solution 
region examined, the net wake skew angle 
appears to be consistent with that predicted by 
momentum. 



(a) Streamlines Released from the Disk 


Centerline 



Comparisons of the computed and experimental 
pressures are shown in the following figures. 
Each figure compares the two prediction 
methods with the experimental data at one 
downstream (X/R) station. The results are plotted 
as pressure coefficients (on an inverted scale) as 
a function of vertical location (Z/R) on the body. 
For reference the vertical location of the rotor 
hub is Z/R = 0.274. 

At the extreme nose of the fuselage (figure 7) 
the AFDD solution matches the experimental 
pressures well. The Mil solution may suffer from 
the coarse panel spacing and the interpolation to 
the measurement location. 



Figure 7: Pressure Comparison at X/R=0.05 

At the second X/R station (figure 8) a better 
match is seen between the two computational 
methods than in figure 7. Both methods show 
less suction over the “shoulder” of the model 
than seen in the experiment. The shoulder is the 
region of the cross-section (shown in figure 3) 
where the radius of curvature is the smallest. 
However, the Mil method does not capture the 
shoulder acceleration at all. This may still be due 
to the coarse panel spacing that does not match 
the curvature of the geometry in this region. 


(b) Vorticity Contours 

Figure 6: Flow Computed from the AFDD 
Solution 
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Figure 8: Pressure Comparison at X/R=0.09 

The trend of both computational methods to 
underpredict the suction at the shoulders of the 
configuration continues at stations up to that 
measured at X/R = 0.31 (figures 8 to 12). 



Figure 9: Pressure Comparison at X/R=0.14 

At more downstream locations on the nose (from 
X/R= 0.14 - 0.35) the Mil method improves in 
predicting the shoulder suction. At X/R = 0.14 
(figure 9) hardly any suction is predicted by the 
Mil method while the AFDD solution shows a 
clear “peak” in the suction similar in trend to the 
experimental data. 



Figure 10: Pressure Comparison at X/R=0.20 

At X/R = 0.20 (figure 10) both methods show a 
peak” in the suction over the shoulder of the 
section. At this station both methods agree well 
with the experiment, except at the bottom of the 
fuselage where the data indicate a deceleration 
of the flow. This deceleration is over-predicted 
by the AFDD method. 



Z/R 


Figure 1 1 : Pressure Comparison at X/R=0.26 

The data from the nose section at X/R = 0.26 
(figure 11) indicates that the flow experiences 
approximately twice the suction at the shoulder 
than predicted by either method. 
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Figure 12: Pressure Comparison at X/R=0.31 

In figure 12 the shoulder suction is matched 
better by the Mil method than by the AFDD 
method, although this may result from the close 
spacing with specific wake filaments at a specific 
azimuth of wake used to produce the Mil average 
influence. 



Figure 13: Pressure Comparison at X/R=0.35 

At the section at X/R = 0.35 (just before the 
nacelle region of the fuselage) both the Mil and 
AFDD methods predict a shoulder suction peak 
between measured data points. At this section 
the AFDD method continues to predict a 
deceleration around the bottom of the section. 
The data at this section indicate a separation, not 
a deceleration at the bottom of the section. In 
fact, at this section the Mil prediction follows the 
trend of the experimental data better than the 
AFDD method at the bottom. The implication 
here is that the Navier-Stokes model does not 


accurately predict the separation that is apparent 
in the data. 



Figure 14: Pressure Comparison at X/R=0.47 

Although no pressure data were obtained on the 
nacelle region of the fuselage, it is clear that the 
two methods disagree in the character of the 
surface pressures over the nacelle at the section 
X/R = 0.47. The prediction of shoulder peak 
suction in the Mil solution for this section suffers 
from coarse panel interpolation, since the 
precise geometry of the nacelle-body junction is 
lost. 



Figure 15: Pressure Comparison at X/R=0.60 

At the X/R = 0.60 section only four experimental 
data points are available for correlation. The 
AFDD solution best fits all four of these points, 
while the Mil solution is a good fit to the lower 
three. Again, similar to the prediction at X/R = 
0.47, the two methods are in significant 
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disagreement on the character of the flow on the 
nacelle of the fuselage. 



Figure 16: Pressure Comparison at X/R=0.88 

Although both predictive methods show that the 
flow at X/Ft = 0.88 is relatively benign, there is a 
discrepancy with the experimental data. The data 
at this station are most affected by the influence 
of the model support. The contribution of the 
support is not included in either solution 
method. 



influence of the fuselage separation seen in the 
experimental data while the Mil method shows a 
stagnation behind the nacelle (approximately 
Z/R of 0.14) and a strong acceleration over the 
shoulder (approximately Z/R of 0.09). 



Figure 18: Pressure Comparison at X/R=1.16 

Behind the nacelle, shown in figure 18, the 
influence of the accelerated rotor wake on the 
surface pressures is demonstrated by the Mil 
method and seen in the experimental pressures. 
Although the magnitudes are not predicted 
accurately, the trend of the Mil prediction is 
good. This overprediction in the magnitude of 
the surface suction may be due to the influence 
of the filaments that are used to model the rotor 
wake and their close influence on the panel 
solution. The AFDD method, in contrast, 
completely misses the acceleration of the flow 
on the shoulder and bottom edge of this 
section. This is, perhaps, due to the numerical 
diffusion of the velocity gradients in the inner 
sheets of the rotor wake by this time-averaged 
numerical method. This diffusion may cause an 
overprediction of the region of separation. 


Figure 17: Pressure Comparison at X/R=1.00 

At locations behind the rotor (X/R > .8) the 
Navier-Stokes method predicts very small and 
relatively uniform values of surface pressures. In 
this region, the Mil method predicts an 
accelerated flow with the lowest pressures on 
the fuselage where there are relatively sharp 
edges that will force flow acceleration. 
Specifically at section X/R = 1 .00, the AFDD 
prediction shows, with surprising accuracy, the 
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Figure 19: Pressure Comparison at X/R=1.35 

Mid-way down the tail region (X/R = 1 .35) the Mil 
method continues to show the strong interaction 
between the wake and the fuselage pressures at 
the top and bottom corners of the fuselage. In 
this region the fuselage transitions from the 
super-ellipse shape with well-defined corners to 
a round section. The experimental data show the 
influence of the wake as an acceleration (or 
suction) over the shoulder. The bottom of the 
fuselage, however, returns to a base pressure 
indicating separation. Here the AFDD method, 
although not accurate in magnitude, matches 
the trend of the experimental pressure. 



Figure 20: Pressure Comparison at X/R=1 .53 

At the final station where experimental data are 
available, X/R = 1.53, only four experimental 
pressure ports are available. Similar to the mid-tail 
station (figure 1 9) the AFDD prediction seems to 
capture the trend of the experimental data. At 
the middle of this section both predictive 


methods have a similar magnitude that is well 
below that of the experimental data. This 
indicates that the mean velocity of the flow at this 
station is underpredicted. 

Concluding Remarks 

A comparison of methods for the prediction of 
fuselage aerodynamics at low speeds has been 
made with one set of experimental data. In 
general, both methods provide insight into this 
complex flow with good general agreement with 
the experimental pressure distribution. From this 
limited comparison, some specific observations 
can be made: 

• The total cost of producing an 
engineering estimate for a rotor-fuselage 
combination is much lower for the model-based 
approach of Mil when compared with the Navier- 
Stokes approach of the AFDD methodology. 

• Although the AFDD method correctly 
predicts the influence of separation behind the 
nacelle region, the apparent separation below 
the fuselage ahead of the nacelle is not captured 
well by this Navier-Stokes method. 

• The choice of breaking the fuselage into 
a specific number of panels will result in some 
inaccuracies in the representation of the airflow. 
This inaccuracy is evident in the comparison of 
the flow on the extreme nose of the fuselage. 

• Over the nacelle, where there is no 
experimental data, the two predictive methods 
show significantly differing flow characteristics. 

• Just aft of the nacelle the acceleration of 
the flow by the wake is seen in the Mil prediction 
of the fuselage pressure, but missed by the 
AFDD method that shows only a region of 
separation. 

NASA Langley Research Center 
Hampton, VA 23681-2199 
April 1, 1998 
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Appendix — Mil Fuselage 
Panelization 

Calculation of an attached flow around a fuselage 
is an important part of the general design of the 
fuselage. For this purpose a method based on a 
theory of potential flows is used. A solution is 
sought as a potential of velocity perturbations. A 
distribution of the potential over the fuselage 
surface is presented as a sum of potentials of a 
double layer (doublets) and a single one 
(sources). 

The fuselage is broken down into a finite number 
of panels and a no penetration boundary 
condition is applied in the center of every panel. 
It is assumed that the distribution of the double 
layer (doublets) and the single one (sources) are 
constant over each of the panels. The density of 
the double layer is determined using a given 
external field of velocities. In calculating the 
surface integrals an assumption is used that the 
fuselage panels are right-angled. Changing from 
the real fuselage to the right-angled panels 
causes distortions at corners of the rectangles. 
However these distortions have little influence 
on the results since the no penetration 
conditions are enforced at centers of the panels. 
This is especially true as the number of the 
panels increases since these inaccuracies will 
grow smaller. As a result, a system of equations 
in unknown potentials of the double layer was 
solved on each of the panels. The solution was 
sought by an iterative method. 

The fuselage breakdown into the panels is a key 
element directly preceding the calculation. The 
breakdown into the panels has to render the 
fuselage shape as exactly as possible on one 
hand and on the other hand to be convenient for 
calculation of velocities on its surface. The 
velocities are derived from an unknown 
distribution of the potential. 

The fuselage breakdown into panels has direct 
influence both on the convergence of the 
iterative process and on the accuracy of the 
results obtained. Programs were developed for 
the fuselage breakdown and for control of the 
resulting panels. These programs locate the 
panel centers, calculate the lengths and 
directions of the normals and tangents to the 
surface, and enable a visual inspection using 
graphical displays. 


A system of coordinates X, Y, Z is assigned to 
the fuselage with the X-axis as longitudinal. The 
fuselage is intersected with planes normal to the 
X-axis. Positions of these planes are selected 
based on specifics of the fuselage geometry. 
The lines of intersection of these planes with the 
fuselage surface create closed contours. 

A relative center may be chosen for each of 
these contours (a center of mass of the cross- 
section may be taken assuming mass to be 
constant over the surface). The contours are 
divided into an equal number of parts. Then 
each of the real contours are projected onto a 
unit circle with its center at the center of mass, 
equal angles onto the real contour of the 
section, and each unit circle is divided into 
subintervals equal in angle. Associated points 
on the real contour can be taken as reference 
points. In cases of the complicated contours, an 
irregular breakdown with unequal angles should 
be used. 

As a result, each of the cross-sections will be 
divided into an equal number of parts. The 
fuselage breakdown into the panels can be 
realized by joining the respective points of the 
adjacent contours. The panels constructed in 
this way are generally not flat. The angles formed 
by their sides are not right. Additional operations 
are carried out to form each of the panels. The 
resulting panel geometry for the ROBIN fuselage 
is shown in figure 4 of this report. 

As a result, radius vectors of the panel centers in 
the adopted system of coordinates, unit normal 
vectors, unit vectors of one of the tangents 
(located in a longitudinal section of the 
fuselage), and lengths of two mutually 
perpendicular directions are computed for each 
of the panels. 

The complexity of the fuselage shape results in a 
disorganized arrangement of the panel centers 
where the potential is to be determined. This 
disorganized arrangement presents difficulties 
for computing the velocity tangent to the surface 
by differentiation of the potential. 

It is convenient to use the fuselage breakdown 
already at our disposal to avoid a double 
approximation of the potential over the whole 
surface. Calculating total velocity at a given point 
requires only derivatives of the potential with 
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respect to length in two different directions. We 
can use the vector of the tangent located in the 
fuselage longitudinal plane as one of the 
directions. A vector located in the plane 
perpendicular to the X-axis and tangent to the 
fuselage contour in the given section is 
conveniently used as the second direction. To 
determine the direction of this tangent, the 
contour at the given point is approximated with a 
circle. 
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